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Abstract. We propose a method to estimate the entropy of thermal pions in 
A+A collisions irrespective of unknown form of freeze-out (isothermal) hyper- 
surface and transverse flows developed. We analyse the average phase-space 
densities and entropies of the thermal pions vs their multiplicities and collision 
energies. The behaviour of these values apparently indicates the deconfinement 
and chiral phase transitions in relativistic nuclear-nucleus collisions. 

Keywords: relativistic heavy ion collisions, phase-space density, entropy. 
PACS: 24.10.Nz, 24.10.Pa, 25.75.-q, 25.75.Gz, 25.75.Ld. 



1. Introduction 

The actual numbers of partonic degrees of freedom released in ultra-relativistic 
nuclear-nucleus collisions are fairly big, up to tens of thousands. As it is commonly 
supposed such a quasi-macroscopic system becomes thermal at the early stage of 
collision process: current phenomenological and pQCD estimates give the times of 
thermalization between 0.6 and 3 fm/c [1]. Then the system expands nearly isoen- 
tropically - this standard point in the hydrodynamic approach to A+A collisions is 
advocated theoretically in the Refs. [2] . Thereby the entropy produced in A+A col- 
lisions carries out the information about very initial thermal stage of these processes. 
It is found recently [3] that the pion phase-space density averaged over momentum 
(totally or at fixed rapidity) and space is also about the constant during the stage 
of the chemically frozen expansion. So, measurements of the entropy and average 
phase-space density (APSD) in thermal hadronic systems at the final (freeze-out) 
stage of A+A collisions gives the possibility to look inside the previous stages such 
as the partonic thermalization and the hadronization (or chemical freeze-out). The 
aim of this work is to estimate the entropy and the APSD of thermal pions in heavy 
ion collisions at the different energies of CERN SPS and BNL RHIC and to analyse 
the physical meaning of the results. 
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2. The extracting entropy of thermal pions 

As well known the entropy of a gas of bosons/fermions has the following covariant 
form 

S =( 2J+1 ) / [-(2*) 3 /ln((27r) 3 /) ± (1 ± (2^) 3 /)ln(l ± (2^) 3 /)] (1) 

where the ± sign corresponds to bosons/fermions with total spin J and a is some 
hypersurface in Minkowski space. The value depends on the distribution function, 
or phase-space density f(x,p) that should be known to make the correspondent 
estimates. It is easy to show, however, that the phase-space density, e.g. of n~ , 
cannot be extracted from two- and many- particle spectra even if the system at the 
kinetic freeze-out is characterized by the locally equilibrium distribution function. 
To study the problem let us write the Wigner function (the analogy of the phase- 
space density for quantum systems) for weekly interacting particles in the mass-shell 
approximation [4]: 

f(x,p) = (2tt)- 3 J d 4 qS(q ■ v)e- iqx {a + (p - (l/2)q)a(p + (l/2)q)) . (2) 

Here the brackets (...) mean the averaging of the product of creation and annihilation 
operators with density matrix refereed to a space-like hypersurfaces where particles 
become or are already nearly free. If the freeze-out is sudden, one uses usually a 
thermal density matrix at the freeze-out hypersurface. The invariant single- and 
double- particle spectra have the form: 

nip) = (a + (p)a(p)) ,n(p 1 ,p 2 ) = n(p 1 )n(p 2 ) + \(a + (p 1 )a(p 2 ))\ 2 , (3) 

and the correlation function is defined as the following 

C(p,q) ^n(p 1 ,p 2 )/n(p 1 )n(p 2 ). (4) 

It is easy to see now that the phase of the two- operator average (a + {p\)a{p 2 )) can- 
not be extracted from the single- and double- particle spectra. It is possible to show 
that the same takes place even when many particle spectra are included into an anal- 
ysis. Therefore one cannot reconstruct the distribution function f(x,p) in a model 
independent way. However, since the f(x,p) is real (note, that for the (locally) 
cquilibritcd quantum systems the Wigner function is positive), one can use Eq. (2) 
to express the f 2 (x,p) integrated over space coordinates just through the squared 
absolute value of the two- operator average, \(a + (p — (l/2)q)a(p + (l/2)g))| , inte- 
grated over q. Then, accounting for the direct links (3), (4) of the later value with 
the single particle spectrum (that is just the f(x,p) integrated over space coordi- 
nates) and correlation function one can extract the phase-space density averaged 
over some hypersurface a at or " after" freeze-out 

iff w Hf( x ^P)) 2 P fld<7 ^ 1 f n ^>) (n i \ 1^3 ( r.\ 

{f{<r,P)) = —fTf = / —Fr{C{p,q) - l)<Tq 5 

J f(x,p)pi 1 da ll (2tt) 3 J p° 
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directly from the experimental data in full accordance with the pioneer Bertsch idea 
[5]. 

At the first sight the extracted value of the APSD does not help to calculate 
the entropy (1), and so, some phenomenological functions that reproduce the ap- 
proximate Gaussian behaviour of the correlation function are utilized usually [6]. 
One should understand, however, that, since we cannot extract the phase of the 
two- operator average, there are an infinite set of distribution functions compatible 
with the observables and, therefore, the entropy calculated will depend on the class 
of the functions we choose which. Here we propose the method how to estimate 
the entropy using just the APSD, do not supposing any concrete expression for the 
phase-space density. 

The method is similar to what was proposed in Ref. [3] to estimate the over- 
population of the phase-space. The idea is based on the standard approach for 
spectra formation [7] that supposes the momentum hadronic spectrum in expand- 
ing locally equilibrated system becomes frozen at some space-time hypersurface with 
uniform temperature and particle number density. Then, within this approximation 
which is probably appropriate in some "boost-invariant" midrapidity interval, the 
phase-space density totally averaged over freeze-out hypersurface and momentum 
excepting for the longitudinal one (rapidity is fixed, e.g., y = 0) will be the same as 
the totally averaged phase-space density in the static homogeneous Bose gas: 

2 

(2.)' </(.,„)U -i§g* (6) 

where 

J eq = (exp(/3(p - M) - (7) 

and P and \x coincide with the inverse of the temperature and chemical potential at 
the freeze-out hypersurface [3] . The essence of the result is that under above condi- 
tions any function of the locally-equilibrium distribution, F(fi, eq ,), integrated with 
the measure p^da^cfip/p contains the common factor, "effective volume" J da^u^, 
that completely absorbs the flows u^(x) and form of hypersurface a(x). This fac- 
torization property has been found first in Ref. [8], multiple used for an analysis of 
particle number ratios (see, e.g., Ref. [9]) and recently of the APSD in Ref. [3]. 

In this work we apply the approach to an analysis of the thermal pion entropy 
per particle, or specific entropy of thermal pions. Using the same approximation of 
the uniform freeze-out temperature and density and Eq. (1) with local equilibrium 
distribution we get the following expression for specific entropy in midrapidity: 



dS/dy _ J d 3 p h/ eg In f eq + (1 + f eq ) ln(l + f eq ) 
dN/dy fd 3 pf eq 



(8) 



In the above ratio due to the factorization property the effective volume is canceled 
and the final expression depends only on the two parameters: temperature and 
chemical potential at the freeze-out. The temperature can be obtained from the 
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fit of the transverse spectra for different particle species and we will use the value 
T = 120 MeV as the typical "average" value for the SPS and RHIC experiments. 
The another parameter, chemical potential, cannot be extracted from the spectra, 
its value could be fairly high even for the thermal pions because of the chemical 
freeze-out and this parameter is crucial for the entropy estimates. We will extract 
it from the APSD analysis following to Ref. [3]. First note that similar to Eq. (5) 
for (f(a,p)) one can express the totally averaged in midrapidity APSD (f(<r, y)) y=0 
by means of the transverse spectrum and correlation function in standard Bertsch- 
Pratt parameterization [5,3]: 

(2«f (/(., y) )y ^ « K V __ J (Q) 

The factor k is introduced to eliminate contribution of short-lived resonances to 
the spectra and interferometry radii [3]. It absorbs also the effect of suppression 
of the correlation function due to the long-lived resonances. Because of the chem- 
ical freeze-out a big part of pions, about a half, are produced by the short-lived 
resonances after thermal freeze-out, that leads to an additional non-equilibrium 
production of the entropy. To eliminate these non-thermal contributions to the 
pion spectra and correlation functions, wc utilize the detailed study of the corre- 
spondent contributions presented in Ref. [3] where within hydrodynamic approach 
it was found that n — 0.65 for SPS and n = 0.7 for RHIC if half of pions is produced 
by resonances at the post freeze-out stage. Then, by means of Eqs. (6), (9) one can 
extract the pion chemical potential. The latter value gives possibility to estimate 
the phase-space overpopulation and the specific entropy of thermal pions as it was 
explained above. 

3. Specific entropy and APSD at SPS and RHIC. 

To evaluate the APSD of direct pions (9) we utilize the following parametrization 
of the 7r~ transverse spectra and interferometry radii. 

For SPS Pb+Pb(Au) 40, 80, and 158 AGeV (^/sjvjv =8.8, 12.3 and 17.3 GeV): 
The transverse spectra are 2 ^m T dm T d y ~ ^ ex P(~ TO T /T e j /), T e ff w 0.169, 
0.179, 0.180 GeV, the midrapidity densities are w 106.1, 140.4, 175.4 (NA49 

Collaboration, [10]). The interferometry radii are Rl = Cl/^Jiiit, Rs = 
Cs,i/ \/l + Cs,2?nT and correspondent numerical parameters are taken from Ref. 
[11] (CERES Collaboration). We use approximation Ro = Rs until minimal mea- 
sured pt momentum, pt = 0.125 GeV, and our analytical approximation of the 
CERES outward interferometry radii data for pt > 0.125 GeV. 
For RHIC Au+Au ^/W = 130 GeV: 

The transverse spectrum is 2ffm ^ Tdy = A (exp(-m T /T e// ) - l)" 1 , T eff w 
0.218 GeV, and ^f- w 249 (STAR Collaboration, [12]). We use here the Bose- 
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Einstein parameterization of the transverse spectra from Ref. [12] since the inte- 
grated rapidity density of negative pions with this fitting function is more close to 
presented recently by the PHENIX Collaboration value ~ 270 at midrapidity 

[13] than what one can get from the exponential parameterization of the transverse 
spectra, dN / y ~ ~ 229 [12]. The phenomenological parameterization of interfcrome- 

try radii, Rl — Cl/^/ttit, Rs — Cs, + Cs,2 + Cs,3m,T , and the correspondent 
numerical parameters are taken from Ref. [14] of the PHENIX Collaboration. As 
one can see from Fig. 3 of Ref. [14], there is some discrepancy in the data on Rs 
radii between the STAR [15] and PHENIX Collaborations [14]. To optimize uncer- 
tainties in the forthcoming estimates , we utilize for the parameter Cs,i the value 
8.75 fm that is the average between "PHENIX data motivated" value, 8.1 fm, and 
"STAR data motivated" value, 9.4 fm, [14]. We use here approximation Ro = Rs- 
The results of our calculations of the APSD according to (9) are used then to 
extract the pion chemical potentials and then the specific entropies basing on Eqs. 
(6) and (8) and supposing the temperature of thermal freeze-out T t h = 120 MeV. 
The values founded are collected in Table 1 and in Figs. 1, 2, 3. We present in 
Table 1 also the freeze-out densities of thermal pions at different SPS and RHIC 
energies. The correspondence between points and energies in Figs. 1, 2, 3 one can 
sec from Table 1. 

TABLE 1 





8.8 GeV 


12.3 GeV 


17.3 GeV 


130 GeV 


dN t l/dy 


106.1/2 


140.4/2 


175.4/2 


249/2 


(2tt) 3 (f(*,y))i: 


0.125 


0.127 


0.137 


0.168 




30.5 MeV 


31.5 MeV 


37.7 MeV 


54 MeV 


dS- h /dy 


213.8 


282.2 


347.3 


471 


dS? h /dy 
dN? h /dy 


4.03 


4.02 


3.96 


3.78 




0.025 fm" 3 


0.025 fm" 3 


0.027 fm" 3 


0.031 fm- 3 



4. Discussion and conclusions 

Here we represent brief analysis and interpretation of the results obtained. As 
was shown in Ref. [3] the pion APSD (f(<J,y)) is approximately preserved during 
isoentropic and chemically frozen expansion unlike fast decrease of the particle n(x) 
and phase-space f(x,p) densities. Thereby the systems reach the kinetic freeze-out 
at some typical particle densities that are similar for different collision energies, 
and, at the same time, their APSD increases with energy since it "memorizes" the 
higher initial densities at higher energies. The rise with energy of the pion APSD at 
AGS is seen from the behaviour of the intcrferometry radii presented in Ref. [11]: 
the correspondent intcrferometry volumes RoRsRl do not change or even drop 
with energy while the pion rapidity densities grow up. From Fig. 1 one can see the 
increase of pion APSD also from SPS to RHIC energies. Only at low SPS energies 
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Fig. 1. The average phase-space 
densities of negative thermal pions 
(27r) 3 (/(<t, y))t h as function of ra- 
pidity density of negative thermal pi- 
ons dN?k /dy at SPS and RHIC en- 
ergies. Table 1 gives correspondence 
between points and collision energies. 

the APSD values are approximately constant. The reason for the latter could be 
the deconfinment happens in that interval of energies. Then one should take into 
account a relatively big gap at that energies between the temperatures of the phase 
transition (hadronization) and chemical freeze-out [16] when the hadronic system 
evolves as chemically equlibrated. While the pion APSD might be fairly high at 
the hadronization stage, it drops during the chemically equilibrated stage of the 
evolution until the system reaches the chemical freeze-out. At higher energies the 
temperature gap between hadronization and chemical freeze-out becomes narrow 
[16] and so the pion APSD becomes frozen and observed value tends to what was 
at the hadronization stage, see points for the SPS highest energy and for the RHIC 
one in Fig. 1. 

The Fig. 2 demonstrates the linear increase of the entropy of thermal negative 
pions with their rapidity density that can be described approximately as dSth/dy = 
28 + 3.6 * dNth/dy. Thereby the specific entropy approaches the massless limit, 
S/N — ► 3.6 at m/T — ► 0, when rapidity density grows, see Fig. 3. The direct 
interpretation might be that the chemically equilibrium pion number is frozen at 
about initial very high temperature that increases with the collision energy. The 
existence of the pions in the quark gluon plasma up to zero binding temperatures 
T w 1.6T C is argued recently [2]. However, the mass of pions in the QGP most 
probably grows with temperature (see, e.g., [17,2]) that turns down the possibility 
of such an explanation. Formally, the decrease of specific entropy of thermal pions 
is conditioned by an increase of the pion chemical potential at kinetic freeze-out. 
The latter happens since the APSD is conserved during the evolution and it grows 
with energy at the chemical freeze-out. Even if the thermodynamic parameters at 
the chemical freeze-out do not change, the APSD grows if the intensity of flows 
increases, - it is determined by the correspondent reduction of the " interfcrometry 
volume" [3]. Since there is the strong correlation between the initial densities in 
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A+A collisions and intensity of the flow developed, the conservation of the APSD 
results in an increasing of observed of the pion chemical potential and an decrease 
of specific entropy with collision energy if the hadronization and chemical freeze-out 
are nearly coincided. 




Fig. 2. The entropy rapidity density 
of negative thermal pions dS^ h /dyas 
function of rapidity density of nega- 
tive thermal pions dN? h /dy at SPS 
and RHIC energies. The solid line 
corresponds to linear approximation 
28 + 3.6 * dN^ /dy. 



dN th /dy 



100 120 140 



d Nth/d y 



Fig. 3. The specific entropy of nega- 



tive thermal pions 



dSJ, /dy 



as func- 



dN t h ' d y 

tion of rapidity density of negative 
thermal pions dNJ/ h /dy at SPS and 
RHIC energies. The solid line corre- 
sponds to specific entropy in massless 
gas. 



Some peculiarities of the energy dependence of specific entropy could be asso- 
ciated also with the deficit of pions at highest SPS and probably RHIC energies 
if one supposes pion chemical equilibrium at chemical freeze-out [9]. The neces- 
sary additional contribution can be associated with decays of a mesons [18]. Due 
to smaller mass of <r at the phase boundary as compare to its vacuum mass, the 
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chiral phase transition is accompanied by an appearance and subsequent decay at 
post hadronization stage of a big number of these particles. Therefore, besides the 
basic increase of the pion APSD with energy at the chemical freeze-out, the above 
mechanism can lead to an increase of the APSD at subsequent evolution and, thus, 
to additional decrease of specific entropy of thermal pions. 
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